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Csl ■ I. ABSTRACT 

o' 

^^ ' Long-term videomicroscopy and computer-aided statistical analysis were used to determine some characteristic 

Qs parameters of in vitro cell motility and proliferation in three established cell lines derived from human glioblastoma 

^\ tumors. Migration and proliferation activities were compared among the three cell lines since these are two features 
"^ 1 of tumor cells that strongly influence the progression of cancer. 

^ ' Cell proliferation in sub-confluent cultures were evaluated by calculating the growth rate of the number of cells and 

' ^ , the distribution of the cell-cycle times in a given microscopic field. In these parameters no significant difference was 
^^" observed among the cell lines regardless to the number of passages. 

f^ , Studies on cell locomotion revealed strong fluctuations in time and exponential distribution of cell velocities. In 

Mh" spite of the fluctuations, both the distribution profile and the average velocity values were reproducibly characteristic 
^ [ to each cell line investigated. 

The results on these dynamical parameters of cell locomotion were compared to pathological data obtained by 
traditional methods. The data indicate that the analysis of cell motility provides more specific information and is 
potentially useful in diagnosis. 



II. INTRODUCTION 

In order to provide appropiate methods for the characterization of the large variety of human brain tumors, repro- 
ducible and quantitative models are needed. Then, one can carry out a characterization of the dynamics of tumor 
cells providing substantial information supplementing the results of such statical methods as histology or molecular 
genetics. For example, in vitro studies on the cell cycle and motility of tumor derived cells could provide relevant 
information for the diagnosis and the estimation of the progression of the cancer. The correlation between tumor 



malignancy and increased migration activity of glioma cells had been demonstrated by in vivo and in vitro exper- 
iments using rat brain xenografts ||l|,|ll| or radial dish assays @j§]. In addition, studies of cell proliferation using 
either proliferation-marker antibodies in tissue sections or flow cytometry of tumor- derived cell cultures have also 
demonstrated that the growth fraction of the specimen often correlates with the malignancy of the tumor |^,|l^ . 

In several cases tumor recurrence is originated by a minor subpopulation of highly migratory cells which were not 
removed surgically with the bulk tumor. Thus, the detection of such subpopulations and the estimation of their 
infiltration depth is essential for postoperational treatments. However, fluctuations in locomotion or proliferation 
generate difficulties in the characterization of the heterogeneity of a tumor specimen; according to early studies a 
strong variation is present in the intermitotic times even within a single-cell derived population lia]. Fluctuations 
dominate the motility of single cells as well |^,|l^,0. The classical methods assaying cell division frequency or cell 
motility cannot distinguish between a fluctuating but homogeneous population and an inhomogeneous population, 
thus their prognostic value is rather limited regarding heterogeneity [10 1. To reveal details about heterogeneity and 
subpopulations in vitro experiments and careful statistical analyses are required. 

Simultaneous and automatized observation of the mitotic and locomotory activity of a large number of cells during 
a long time interval can provide a direct way to overcome such difficulties. 

We describe here a novel long-term videomicroscopic system together with a computer-aided statistical analysis 
to quantify the in vitro behavior of cells. With this method we compared three established glioblastoma cell lines 
derived from three patients. Glioblastoma is a common and highly malignant (WHO Grade IV.) tumor that develops 
in late adult life and generally located in the cerebral hemispheres. These tumors contain motile and invading cells, 
frequently causing rapid recurrence after surgical resection |l5[ . In some cases the migration of these cells lead to the 
progress of the disease even without the formation of notable mass effect [|9| . 

While our studies revealed a rather similar proliferation activity in all investigated cell lines, a reproducible and 
significant difference was found in the locomotory activity of the cells. Due to the large number of individually tracked 
cells the question on the homogeneity of cell motion and proliferation of the cell populations could be also addressed. 
By appropriate automatization the complete analysis can be completed within 4 or 5 days after processing the surgical 
specimen to one-cell suspension, providing a new tool for diagnostical purposes. 

III. MATERIALS AND METHODS 

A. Diagnosis of human brain tumors 

The surgical specimens were obtained during craniotomies for resection of 3 hemispherical primary brain tumor. The 
bigger part of the specimen was evaluated by routine histopathology, including staining with hematoxylin and eosin 
and immunohistochemical staining for GFAP expression. Specimens were graded according to the WHO classifications 
|l2| at the Histopathological Department of National Institute of Neurosurgery. 

B. Establishment of human brain tumor cell lines 

Part of the specimen was used for establishing the cell lines. The samples were washed by Minimal Essential 
Medium (Sigma) containing gentamycin (Chinoin) and Fungizone (Gibco). After removing the vessels, the samples 
were minced and then triturated with a Pasteur-pipette. The primary tissue was seeded into Leighton culture flask 
on Bellco-slide or into Steriline culture flask in Dulbecco's MEM (Sigma) with 20% FCS (Gibco). According to our 
procedure the cultures are not one-cell derived. The passages were performed with 0.25% trypsin and cells transferred 
to Greiner tissue culture flasks. Cells were maintained in DMEM (Sigma) with 10% FCS, 40 /^g/ml gentamycin and 
5 Aig/ml Fungizone, at 37°C in a humidified 5% CO2 atmosphere. Cells were stored in 10% DMSO (Sigma) DMEM 
solution in liquid nitrogen for later use. 

C. Immunocytochemical staining for GFAP in the cultures 

The cultured cells' glial origin was indicated by the presence of glial fibrillary acidic protein (GFAP). Cells cultured 
on slides were fixed with Zamboni-solution for 60 minutes at room temperature. Then incubated as follows: 3 times 10 
minutes in PBS (Phosphate Buffered Saline), 20 minutes 15% human albumin in PBS (Human), 30 minutes with mouse 
anti-GFAP antibodies (High Performance, Biogenex) in a humidified chamber, 20 minutes biotin-conjugated anti- 
mouse immunoglobulin (Multilink Stravigen, Biogenex), 20 minutes peroxidase-conjugated streptavidin (Biogenex). 



Slides were developed with DAB (2,5 mg DAB+ 5 ml PBS+ 50/il 1% H2O2) under microscopic control. Controls 
were performed without the primary antibodies. 

D. Long-term videomicroscopy 

10^ cells, counted in haemocytometer, were plated in 35 mm TC-dishes (Greiner) with 2 ml medium (DMEM 
(Sigma) with 10% FCS (Gibco)). Cell cultures were kept in a mini-incubator - providing 37°C in a humidified 5% 
CO2 atmosphere - attached to the powered stage of an inverted phase-contrast microscope (Zeiss Televal-l)(Fig 1). 
Images of 3-6 neighbouring microscopic fields were taken in every 5 minutes during a 3-day long period, with 10 x, 
20 X or 40 X objectives using a CCD camera (JVC KY-F30B) connected directly to the frame grabber card (Matrox 
Meteor, Matrox Electronic Systems LTD, Canada) of a PC (running under LINUX operating system). 

E. Cell positions and trajectories 

To determine the position {xl) of the individual cells, the geometrical center of each cell was tracked manually in 
every 4th image (i.e., in every 20 minutes in real time). The difference in the location of a given cell was considered 
as a migratory segment. Trajectories were constructed from these segments. 

F. Duplication time 

Using the database of cell positions, the total number of cells in a given microscopic field (n) was determined. The 
growth of cell number was found to be approximately exponential in each investigated culture, i.e., n ~ e"*. Based 
on a least square-fit of a, the duplication time, r, was calculated as r = (In2)/Q!. 

G. Cell-cycle time and non-proliferating cells 

All cell divisions during the observation period were identified resulting a data base containing the corresponding 
mother and daughter cells. The cell cycle length (tq) of a given cell was then determined as the time elapsed between 
two consecutive cell divisons. 

We denote by v the rate at which non-proliferating cells appear in the culture. It was calculated in two steps. 
First we determined N, the number of proliferating cells by the criteria that both their birth and subsequent mitosis 
was recorded and their birth occured more than 30 hours before the end of the recording. Then N^ , the number of 
suspected non-dividing cells was calculated. These cells were selected using the following criteria: (i) lack of observed 
mitosis during the entire recording period, and (ii) being tracked for at least 30 hours following their birth. Then, i> 
was calculated as 1/ = N^,/{N^, + N). 

Note, that there are two sources of systematic errors in this procedure. On one hand, when calculating A^*, cells 
may be included which do divide, but their cell cycle length is longer than the time period during the culture was 
observed. On the other hand, we did not include in A^^ cells which do not divide, but we could not track them for 
longer than 30 hours because they left the field of observation. Taking into account that during a 30-hour period 
approximately 10% of the cells migrated out of the observed area, and (based on the empirical distribution of the 
cell-cycle times) « 20% of the dividing cells have longer cell cycle than 30h, we estimate the relative error of v to be 
10%. 

Assuming that after each division a J^th portion of daughter cells become unable to further proliferate, the size of 
the cell population (n) after the kth complete set of cell division - denoted as n(A;ro) - is given by 

log n{kTo) fn k log 2{1 — u) + const. (1) 

Thus based on this simple calculation, the relation 

Tolog2 = Tlog2(l-i.) (2) 

is expected to hold among the independently measurable parameters t,to and v. 



H. Cellular velocities 

The velocity, Vk{t), of a given cell k was calculated as the translocation of its geometrical center during a given 
time interval, At, i.e.. 



Vkit) 



Xk{t + At) - xk{t) 



At 



(3) 



At was chosen to be Ih, since during this time interval the displacement of the cells were typically larger than the 
cell size, thus the relative error of Vk{t) resulting from the tracking procedure was reduced. Nevertheless, during one 
hour (as our data show) the fluctuations in cell motility were not yet averaged out. 

To characterize such a fluctuating quantity as the cell locomotion activity we determined the F(v) cumulative 
distribution function. F{v) gives the probability of the event that the velocity of a randomly selected cell at a given 
time is less than v. Since the relative error of the small cell displacements is high due to the manual tracking procedure, 
velocities less than 5 /im/h were discarded in the analyses. The F{v) empirical distribution function was fitted by an 
exponential distribution [Q as 

Fiv) = 1 - expi-v/vo), (4) 

where Vo is a fitting parameter being equal to the mean velocity if (0) is exact. 

As another way to characterize the motility of the cells, the empirical average velocity of the culture, V, was 
calculated as 

y ^ J2kJ2^Mt^) ,g. 

J2kJ2^^kit^)' 

where Sk{ti) is 1, if the fcth cell is in the visual field in the ti moment, otherwise 0. 

IV. RESULTS 
A. Histological description of tumors and cell lines 

The main features of the cell lines are given in Table 1 . The routine histology of the tumor specimens concerning 
WHO classification revealed glioblastoma. All three sections showed hypercellularity, pleomorphism, vascular prolif- 
eration and necrosis. The HA and HC specimens (detailed in Table 1) were classified as glioblastoma multiforme. In 
the HB specimen some areas were dominated by multinucleated giant cells and the invasion of lymphoid cells was 
observed. Thus HB was classified as giant cell variant of glioblastoma. 

The cultures were morphologically characterized according to Bigner . Usually the HA-derived cells had two or 
three processes and the formation of filopodia was quite frequent. The HB culture's cells were mostly elongated, had 
more processes and the formation of filopodia and lamellipodia with ruffling edges was exceedingly intensive. Both 
morphology was described as fibroblastic. In contrast, even the subconfluent HC cultures showed typical epithelial 
morphology with the general absence of cellular processes. The cells formed preferentially a monolayer, though cells 
could crawl over each other. During cultivation multinucleated monstrocells were observed in all the three lines. On 
the time-lapse records cell divisions into three daughter cells were also observed occasionally. In most of these cases 
the resulting daughter cells seemed to perform normal mitoses later on. 

B. Cell proliferation 

The number of cells in a given microscopic field grew exponentially after an initial lag phase of 15-20 hours (Fig 2). 
Based on these data the cell duplication times were determined (Tab 2). No significant differences were found in the 
duplication time of the various cell lines: each culture could be characterized by r= 38±4 hour. 

The cell number in a given microscopic field changed because of cell divisions and the migration of the cells. 
According to our observations the net current of migration was approximately zero, so the increase of cell number 
was basically determined by the ratio of dividing cells and the duration of their cell-cycle. 

Continuous observation of the cultures allowed to determine the cell-cycle time (tq) of individual cells as well. 
The empirical distribution functions of the cell cycle lengths were similar in all the three cell lines investigated (Fig 



3). In fact, based on this amount of data the distributions could not be distinguished based on a Wilcoxon test 
(with p<0.05). The average ceh-cycle time was found to be 25.6±6.2 hours. The difference between the average 
cell-cycle time and the duplication time of the cell population can be explained by the presence of non-dividing cells. 
The proportion of these cells, v, was found to be 32 ± 5% (Table 2). According to relation (||) this proportion was 
determined by average cell-cycle time (tq) and duplication time (r) as 34%. Thus, relation (ph holds within error for 
the rate of arising of non-proliferating cells, average cell-cycle time (tq) and duplication time (r). 

To investigate the role of hereditary factors in the determination of cell-cycle time the correlation between parent 
and daughter cells was studied (Fig 4). The scattering of data indicates uncorrelated variables, it means that the 
duration of cell-cycle is not influenced by the parent cell's cell-cycle time. 

C. Cell locomotion 

The dynamics of cell shape can be qualitatively observed and characterized using the long term time-lapse 
records(Fig 5). The videomicroscopic images showed the rich dynamics of the formation of lamellipodia and filopodia. 
The intensity of ruffling edge formation correlated with the motility of cells ||l^. As expected, the epithelial HC 
cells displayed the lowest motility and the HB cells were more motile than the other fibroblastic HA cell line. 

The analysis of the pathways of individual cells revealed that in subconfluent cultures the direction of cell movement 
is random and disordered (Fig 6a). However, at high cell densities the HB cell line displayed a rather ordered collective 
migration (Fig 6b). 

In all of the three cell lines the cell velocities were highly fluctuating as demonstrated in Fig 7. 

The F{v) cumulative distribution functions are shown in Fig 8. Both the distribution functions and the average 
velocities of the culture {V and uo) were different among the three cell lines but within error reproducible for each 
individual cell line (Table 3). 

D. Population homogeneity 

To understand the relation between velocity distribution function and the homogeneity of locomotory activity of 
the cell population further statistical analysis is required. The velocity distributions described by (|j) were found 
for relatively large (« 100 cells) populations. On the level of individual cells the exponential behaviour (0) can be 
interpreted in two ways: (i) The culture is inhomogeneous, i.e., slower and faster cells can be distinguished on the bases 
of well preserved phenotypic properties. In this case the exponential F(v) distribution can reflect the ratio of the slow 
and fast cells in the culture, while the velocity fluctuations of the individual cells can show an arbitrary distribution, 
(ii) If the culture is homogeneous, then almost all cells exhibit the same distribution of velocity fluctuations, i.e., 
Fi{v) « F{v) holds for each cell i, where Fi{v) denotes the velocity distribution function of the cell i. In this scenario 
the time-averaged velocity {vi{t)) of each cell would be the same if we could calculate the time averages over an 
infinitely long time. Since the time averages are calculated over a finite time T only, for the distribution of the 
average velocities we can expect a Gamma distribution with a parameter s = T/tg, where the correlation time of the 
process is denoted by io Q]- 

We used the distribution function of the average velocity of individual cells (averaged over the entire observation 
period) to characterize the inhomogeneity within the cell lines. In all cases the average velocity of the cells showed a 
gamma distribution indicating the homogcnity of the cell population (Fig. 9.). However, if applying the same kind 
of analysis on a cell position data base in which cells from different cell lines were mixed, a significant deviation from 
the gamma distribution can be observed (shown in the insert of Fig. 9.). 

V. DISCUSSION 

The comparison of classical morphological data and the dynamical properties in the three glioblastoma cell lines gave 
some unexpected results. In spite of the different origin and different passage-levels of the investigated lines significant 
alterations in the in vitro proliferating potential were not revealed. Both the cell-cycle time and the percentage of 
dividing cells under the same culture conditions were found equal in the three different glioblastoma lines. The cell- 
cycle time measured in our cultures did not differ from the doubling time of normal primary fibroblasts from human 
skin p] . The variability of the cell-cycle time of individual cells seems to be independent from inheritable factors and 
might be best explained by the stochastic activity of the cells. 



The experimental verification of Eq. (|^) suggests that the non-dividing cells did not belong to an initial subpop- 
ulation in the culture but arose continuously during proliferation. The considerable proportion of these cells may 
indicate that defective mitoses are abundant in the cell culture. The cause of this phenomenon can be the genetical 
instability of tumor cells which has been recently demonstrated Q . 

Although the in vitro cell motility is a highly stochastic process regarding the direction of the movement or the 
fluctuation of the velocity of an individual cell, the locomotory activity of the cell lines have some quantitative 
properties that can be interdependent with the morphology of the cultured cells. The cell lines with fibroblastic 
morphology performed higher locomotory activity. Reproducible and significant difference was found in the locomotion 
of the cells among the three investigated glioblastoma cell lines. 

The HB case can illustrate that classical studies are not always satisfying in the forecast of the disease course. In 
spite of the diagnosis as a giant cell variant, which is commonly known as a less malignant type of glioblastoma, and 
the powerful presence of immune response, the tumor was fatally renewed within 4 month. It can be supposed that 
the fast and strong recurrence was due to the highly motile cells that invaded the surrounding tissue and were not 
removed surgically. 

Our data were obtained on established cell lines after multiple passages, in this way these cultures may not reflect 
entirely the properties of the cells of the investigated tumors. With similar preparations the evaluation can be 
performed on primary cultures with no passages as well. Investigations of in vitro dynamic behavior of primary 
cell cultures from surgical specimens can complement classical diagnostic methods. Although during culturing an 
inevitable selection takes place, some of the initial features or even heterogenity of the tumor tissue may be reflected 
in the cell cultures. If such analysis can be carried out within a short period of time, it can help to choose the most 
appropriate postoperative therapy in each individual case. 
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Figure legends: 

FIG. 1. Schematic representation of the long-term videomicroscopy system. The incubator, placed on the inverted micro- 
scope's stage, holds the TC-dish. The microscope stage and the adequate cell culture conditions (temperature, atmosphere) 
are controlled by the computer. The CCD camera's images are transmitted to the computer, where the digitalized images can 
be processed and analysed or recorded on videotape. 

FIG. 2. The growth of the normalized cell number in a given visual field. The y-axis is logarithmic, and the actual cell 
numbers were normalized by the A^o coefficient of the exponential fit N ^ Noe"^. The growth rate is similar in all the cell 
cultures investigated. The solid line shows the fitted exponential growth with the average duplication time (r= 38±4). 

FIG. 3. The distribution of the cell-cycle time (ro). The distribution functions obtained from the three cell lines could not 
be distinguished statistically. Thus the proliferating cells in any of the investigated cultures divided with the similar cell-cycle 
time. 

FIG. 4. The daughter cells' cell-cycle time (ro(d)) versus the parent cells' cell-cycle time (ro(p)). The scattering of the data 
points means that the cell-cycle times in the subsequent generations do not correlate. 

FIG. 5. Snapshots from the recordings of the various cell lines. Despite the same diagnosis, morphological differences can 
be observed. The approximately epithelial morphology of HC (c) cells can be compared to the fibroblastic characteristics of 
HA (a) and HB (b) cells. Field of view is 740 //m x 560 pm. 

FIG. 6. Cellular trajectories in the HB3 measurement, (a) During the first 24 hours the cells perform a persistent random 
walk, (b) In the second 24 hours, at a higher cell density, the cells show a spatially ordered migration with increased persistence 
length. Field of view is 740 fim * 560 /^m. 

FIG. 7. Cell velocities fiuctuate in time. In the figure 3 representative cells' velocities are plotted vs time for each cell line. 
For better visualisaton the consecutive curves are shifted vertically by 50/im/h. 

FIG. 8. The distribution function G{v) = 1 — F{v) of the cell velocities in the investigated cultures on linear-logarithmic 
plot. The linearity of the curves indicate that the cell velocities can be described by an exponential distribution. The graph 
also shows that G{v) function rcproducibly distinguish the three different cell lines. 

FIG. 9. Distribution functions of the average velocities of individual cells. The average was calculated over the entire 
measurement period and the functions are normalized by the average velocities of the cultures. The curves are well fitted 
by the calculated gamma distribution suggesting a phenotypically homogeneous population with respect to the locomotory 
activity. In the inserted graph the mixed curve - data from two different cell lines - shows a significant difference from the 
calculated gamma function. 



TABLE I. The most important data of the cell cultures analyzed 



measurement 


HAl 


HA2 


HB3 


HB4 


HCO 


HC2 


date of 
operation 


1997 
March 


1997 
March 


1996 
October 


1996 
October 


1995 
February 


1995 
February 


age and sex 
of patient 


58 y 
f 


58 y 
f 


51 y 
m 


51 y 
m 


61 y 

f 


61 y 

f 


passage 


7 


8 


45 


47 


103 


107 


area of 

observed 

field [mm] 


1.14 


0.41 


0.41 


1.14 


1.14 


0.41 


initial cell 
rmmber 


55 


27 


61 


110 


93 


33 


initial cell 
density [mm]^^ 


48 


65 


148 


96 


82 


80 





TABL 


il II. The data on proliferation in the six measurements 






measurement 


HAl 


HA2 


HB3 


HB4 


HCO 


HC2 


average 


average cell- cycle 
time and its standard 
deviation [h] 


24.2 ± 6.6 


27.9 ± 5.8 


n.a. 


24.7 ± 5.9 


n. a. 


25.9 ± 5.9 


25.6 ± 6.2 


duplication time[h] 


39,6 


33,8 


42,0 


40,77 


43.1 


33,8 


38 ±4 


rate of non- 
proliferating cells 


39% 


n.a. 


n.a. 


24% 


n.a. 


30% 


32 ± 5% 



TABLE III. The average velocities in the six measurements 



measurement 


HAl 


HA2 


HB3 


HB4 


HCO 


HC2 


V bm/h] 


11,01 


9,19 


15,6 


18,6 


4,2 


6,4 


vo ^m/h 


11,6 


11,1 


15,6 


17,2 


6,5 


7,0 


error [/im/h] 


1,2 


1,1 


1,6 


1,7 


0,7 


0,7 



4 - 



3 - 



N/Ni 







1 1 
- 'hal' ■ ■ ■ - 


1 1 


1 




'ha2' 






J 


'hb3' ■ ■ ■ ■ 






^y/ 


'hb4' 




J- 


-^ J 


- 'hcO' ■ ■ ■ 




y 


n> 


'hc2' 




J f 


O^ 


'fitted' 




^ 




^__. 


j^^y'^f'-'"- 


.m'x 




7C^ 


1 1 


1 









10 



20 



30 

t\h\ 



40 



50 



60 



10 




To 



fl 



11 




20 / X 25 
Mp) [h] 



12 






1 


i 


J 




J«-l-' . 1 M 




'•^ 




13 





14 



V 

[/im/h] 




15 



0.1 



G{v) 



0.01 



0.001 



v\ 


1 1 

'hal' 

'ha2' ■ ■ ■ - 

'hb3' 

'hb4' ■ ■ ■ ■ 

'hcO' 

'hc2' ■ ■ ■ 


Vn 


V,^ 







20 40 60 

V [fim/h] 



80 



100 



16 



1 



0.8 - 



FivlvQ 



0.6 - 



0.4- 



0.2 - 








'hal' — 

'ha2' ■ ■ 

'lib3' — 

'hb4' ■ • 

'hcO' — 

'hc2' ■ 

'gamma' O 



1/1.5 2 



2.5 



17 



